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Quantum mechanical systems of identical particles provide an elegant situation in which 
one can apply the ideas of path integrals in multiply connected configuration spaces 
The reason is clear if one permits oneself to ignore coincident points in the configuration 
space, i.e., configurations where two or more of the particles are at the same position ||. In 
that case, non-contractible closed loops in configuration space can easily be constructed. For 
example, a continuous evolution of the system involving an interchange of particles forms a 
closed path due to the identical nature of the particles, yet the path is non-contractible due 
to the interchange. 

According to the general rules of path integrals in multiply connected configuration 
spacefj one can divide the space of paths into homotopy classes, and rewrite the path 
integral as a sum of sub-integrals, each of which is a path integral over one such class: 

f Vx(t)e iS = J2 I Vx(t)e lS } (1) 

where C n denotes the homotopy class labelled by n. Since the main requirement of the path 
integral is that it satisfy the Schroedinger equation, and since each sub-integral in (|1|) is itself 
a solution, any linear combination of the sub-integrals will still satisfy the Schroedinger 
equation. For consistency of the path integral, we must impose one requirement on the 
coefficients: that they form a commutative representation of the first homotopy group of 
the configuration space. 

A very instructive example of the above ideas, due to its simplicity, is provided by a free 
particle constrained to move on a circle There, the paths are classified according to 

the net number of windings around the circle. The topology involved can be regarded from 
the point of view of the simply-connected covering space of the configuration space, which 
is simply a line. Each point on the circle is the image of an infinity of points on the line, 
and different homotopy classes on the circle correspond to classes of paths on the line with 
different pre-images of the final point under the mapping between the circle and line. The 
sub-integral corresponding to any homotopy class on the circle, viewed in the covering space, 
is an ordinary path integral between the initial point and the corresponding final point on 



1 For an excellent discussion, see Ref. [Q]. 
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the line. The path integral for the circle, then, becomes an infinite sum of free-particle path 
integrals in one dimension. 

One can, as described above, add relative phases exp —ina, where n denotes the number 
of times a path winds around the circle, to the sub-integrals; such liberty, in fact, corresponds 
to threading a magnetic flux $ = a/e through the circle 

Applying these ideas to systems of identical particles, the above reasoning was used to 
demonstrate that only Bose and Fermi statistics are possible in three space dimensions |§. 
The reason is essentially that the first homotopy group of the configuration space of N 
identical particles in 3 dimensions is the permutation group Sn, which has only two com- 
mutative representations: the trivial representation and that representation which assigns 
±1 to a path according to whether it corresponds to an even or odd permutation of the 
particles. The former describes Bose statistics, the latter, Fermi statistics. 

The work of Ref. 0] can in fact be turned around slightly: using a path integral formu- 
lation, a system of identical spinless fermions could be described by bosons, by dividing the 
path integral into its homotopy classes and adding a factor —1 to all paths corresponding to 
odd permutations of the particles. Although this is true in principle, the division of the path 
integral into homotopy classes appears to be rather difficult, and to our knowledge such an 
approach has not been put into practice. 

More recently, it was realized by completely independent means that fractional statistics 
(anyons) is allowed in two space dimensions A subsequent re-examination of the 

ideas of Ref. [§] in two space dimensions || revealed in an elegant way that indeed the 
possibility of fractional statistics might have been anticipated via path integrals: the relevant 
homotopy group is a braid group (rather than a permutation group), and the commutative 
representations of such groups allow for any relative phase between different homotopy 
classes. 

Anyons are normally described by adding fictitious ("statistical") charges and magnetic 
fluxes to "normal" particles (bosons or fermions) in such a way that no new forces are intro- 
duced but so that windings of each particle around the others are counted in an appropriate 



2 For a comprehensive review including many applications, see Ref. [||. 
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way. As in three dimensions, a path integral approach provides, in principle, an alternative 
means of altering the statistics. One can imagine dividing the path integral for identical 
bosons into homotopy classes, and adding relative phases between different classes to change 
the way different windings interfere with one another in a way appropriate for the desired 
statistics. 

Once again, putting this program into practice is extremely difficult. The problem is that 
one needs a set of coordinates which retain a memory of the relative windings of the particles. 
(For instance, Cartesian coordinates are clearly inappropriate since the coordinates of a given 
configuration of particles are unique; the winding of particles around one another is lost in the 
Cartesian description of the motion.) For the case of two particles the resolution is obvious: 
if, after disposing of trivial center-of-mass motion, the relative motion is described in polar 
coordinates, the angular coordinate measures the winding of the two particles around one 
another. The path integral can thus be divided into homotopy sectors where the final point 
in configuration space is described by final angles differing by a multiple of it. 

Indeed, a program along these lines has been put into practice in the context of a path 



integral description of the Aharonov-Bohm effect |I0| , |TT| , which has obvious and well-known 
similarities with anyons M, and in the context of the statistical mechanics of anyons JT2 



However, these discussions of the Aharonov-Bohm effect are in a sense "less topological" 
than, for instance, the treatment outlined above of the path integral for a particle on a 
circle Essentially, a (non-topological) path integral in polar coordinates is written 

down for a particle moving in the plane [l3|jljj], and a topological constraint on the path 
integral is added after the fact. The intermediate angular positions are integrated only over 
the interval [0, 2tt], in contrast with the general approach discussed in Refs. [0,0. 

In this paper, we present an alternative approach to the problem of the relative motion 
of two anyons, which has a certain aesthetic advantage over that described above, in that 
it is more faithful to the topology of the situation. The approach consists in regarding the 
motion from the point of view of the simply-connected covering space for the plane with the 
origin excised, which can be visualized either as a sort of spiral staircase or as a half plane, 
and performing the path integral on this space, allowing the angle of intermediate steps in 
a path to take on any real values in a manner exactly analogous to Refs. [TLB], rather than 
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restricting them to the range [0, 2tt]. Our approach is similar in spirit to that of Khandekar, 



Bhagwat and Wiegel [HJ , who discussed the problem of the motion of a particle in a plane 
with the origin removed from an equally topological point of view, although the details of 
our calculation differ from theirs. 

We begin with the Lagrangian describing the relative motion of two bosons of unit mass 
moving in the plane with the origin removed, in polar coordinates: 

L = \(r* + rW). (2) 

Rather than adding an explicit interaction term to (§) to alter the statistics of the particles, 
we will introduce relative phases between the different homotopy classes of the path integral, 
as described above. 

The propagator for motion from a point r' = (r', 6') to r" = (r", 9") in a time interval T 
is the sum over all paths from the initial to the final point: 

Kir"; r'; T) = f Vrit) e iS . (3) 
Jr' 

If we write the paths in polar coordinates, (r(£), #(£)), then (r(0),#(0)) = (r',8'), and 
(r(T),9(T)) = (r",9" + rnr) = r", reflecting the fact that changing the final angle by an 
integer multiple of it describes the same final state: the paths between a given initial and 
final configuration can be divided into homotopy classes labelled by n. 
We can therefore rewrite the path integral in the following way: 

ir(rV^T) = W%r(t) (4) 

n Jr n 

We are now, according to Ref. |J, free to add phases between the sub-integrals in (Q): 

K^K ia) = Y / e~ ina K n - (5) 

n 

The phases change the interference between different homotopy sectors, reflecting altered 
statistics: if a = n (or any odd multiple of 7r), the amplitudes of paths whose n's differ by 
one subtract rather than adding, as appropriate for Fermi statistics, while if a is any even 
multiple of tt, the additional phases in (|5[) have no effect, and the particles are bosons. For 
arbitrary a, the particles are anyons. 



We see that anyons can be described by perfectly conventional path integrals within each 
homotopy sector (i.e., without introducing the statistical interaction), if the total propagator 
is computed with appropriate phases between the different sectors. 

Let us compute, then, a sub-integral, K n . The computation is performed in polar co- 
ordinates, which introduces certain subtleties. In particular, it was shown in Ref. [[13] that 
the most naive guess for the path integral must be modified, due to a path integral version 



of operator ordering ambiguities. In Ref. it was shown that a path integral in curved 



space requires the inclusion of an additional "effective potential" term to the Hamiltonian; 
applying that result to the case at hand (where polar coordinates introduce a nontrivial 
metric, so the formalism of Ref. [16| applies directly), this effective potential is —1/(87 
and the sub-integral becomes 



ir 2 ) 



i /■ N TT 1 ,. J „ J ^„..„..#-A''i-^i) 2 +ft-m " 



where the limit iV — > oo will be understood throughout, 6q = 9', 9^ = 9" + mr, and fj is the 
average position for the interval (J — l,j), which we take to be the geometric mean y/rjZifj. 

We can simplify the angular integrals considerably by the following trick. First, we insert 
a factor 1 = / (19^5(9^ — (9" + nir)). Second, we rewrite the argument of the 5-function in 
terms of the angular differences A9j = 9j — 9j-x, it becomes A9i + . . . + A9n — 59 n , where 
S9 n = 9" + mi — 9' is the total angular change for the propagator for the nth sector. Third, 
we rewrite the Dirac 5-function as an exponential: 

5 A% - 89^ = J ^ exp zA (E A9j - 59^ . (7) 

The N angular integrals can now be rewritten in terms of integrals over the N angular 
differences A9j; the integrals so obtained are independant Gaussians and we obtain 

1 f dX^ase n t^u ^ x A„ < £ (i^-r^f (A 2 - l/4)e\ 



= / j n w n(^)- 2 \ e 



fo) 2 



(8) 



The integrand can be usefully rewritten in terms of modified Bessel functions, using the 
asymptotic form for these; one finds 

i(r' 2 +r" 2 )/2e . J\ .N-l N N-l 

K n = W j ge"*^ J n fadr,) n I W ((fjf/ie) exp i £ rf/e. (9) 



The radial integrals can now be performed ||14]]TTH , yielding (taking the limit N — > oo 



r " A -i\&e n 



Kn = ^ J ^J W (-i^/n (10) 

our final result for the sub-integral. 

Although we have derived this result for a system of two identical bosons, it applies 
equally well to the problem of a single particle moving in a punctured plane (which is also, 
of course, described by the Lagrangian (|2|), if we write 59 n = 9" + 2rm — 9', since for a single 
particle we have periodicity of 2tt rather than tt. With this change, our result fllQ|) is then 



in complete agreement with that of Ref. [15 



A simple check of this result is to sum the sub-integrals (taking 59 n = 9" + 2n-n — 9') 
for all the topological sectors; one must obtain the propagator for a free particle. Instead, 
we will perform this sum with relative phases added in, thus obtaining the propagator for 
a charged particle in the presence of a flux tube; the free particle is obtained in the limit 
where the flux is zero. The propagator (for e<p = a) is 

Kia) =T,e~ ina K n = — J ^^ e -^ +a ^ 2 ™e-^ e ''- e 'H w (-ir'r"/T). (11) 

Shifting A — > \ — a/2n and using the Poisson summation formula J2 n e~ lX2nn = J2 m cx(X+m), 

i(r' 2 +r" 2 )/2T 

= £ c'^*^-^ W/^iC-i^/T). (12) 

L'Kll m 

This is our final result for the propagator in the presence of a flux tube. It differs from 
the result of, e.g., Ref. ||11||P| by an overall phase exp i(a/27i)(9" — 9'). This phase is in fact 



simply a choice of gauge for the problem |T0[, and is related to whether or not one chooses 
to describe the particle by a single- valued wave function (as in |[L1|| ) or not (as in this work). 
Alternatively, our method of describing the interaction of the particle with the flux tube 
takes into account the additional phase for paths which wind around the solenoid different 
numbers of times, but it does not take into account the phase due to the change in angle 
between the initial and final points; if we wanted, we could have taken this into account by 



3 Note that Ref. Ol uses a definition of a which differs by a factor 2tt from ours. 
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adding an overall phase exp —ia(9" — 9')/2tt to the propagator, in which case (|12|) would be 
in complete agreement with Ref. | |11|| . 

We can take the limit a — > in the sum is easily evaluated in terms of the 

generating function for modified Bessel functions and one finds that the free propagator 
is indeed obtained. 

For the application of these ideas to the problem of the relative motion of two anyons, 
the change is trivial: different homotopy sectors now differ in their final angle by multiples 
of 7r (rather than 2ir), and the change in angle for the nth sector is now 9 n = 9" + rm — 9'. 
The propagator becomes 

p i(r' 2 +r" 2 )/2T 

= - £ e^ 2m+ ^( e "- e ')/| 2m+a/ .|(-^V7T). (13) 

ml m 

As expected, one can easily verify that when a = or 7r (corresponding to bosons or 
fermions) the propagator is 

2mT V I v ' 

Unfortunately, applying these ideas to systems whose configuration space is more compli- 
cated {e.g., three or more anyons), while in principle straightforward, appears exceedingly 
difficult because of the difficulty in finding a suitable set of coordinates which retains a 
memory of the windings of particles around one another. While such coordinates have been 
found for the case of three anyons, and in fact an expression analogous to (ffTl) can be found, 
in those coordinates the action is sufficiently horrible that the path integral does not seem 
feasible. 

In summary, we have given a derivation of the propagator for a particle in the plane 
with origin removed which is in agreement with previous work, but which is more faithful 
to the topology of the configuration space. We have added permitted phases to different 
sectors of the path integral in order to describe the Aharonov-Bohm effect and a system of 
two anyons. 

This work was supported by the Natural Science and Engineering Research Council of 
Canada. 
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